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Abstract; The Nekrasov-Shatashvili limit for the low-energy behavior oi Af = 2 and Af = 2* 
supersymmetric SU{2) gauge theories is encoded in the spectrum of the Mathieu and Lame equations, 
respectively. This correspondence is usually expressed via an all-orders Bohr-Sommerfeld relation, but 
this neglects non-perturbative effects, the nature of which is very different in the electric, magnetic 
and dyonic regions. In the gauge theory dyonic region the spectral expansions are divergent, and 
indeed are not Borel-summable, so they are more properly described by resurgent trans-series in 
which perturbative and non-perturbative effects are deeply entwined. In the gauge theory electric 
region the spectral expansions are convergent, but nevertheless there are non-perturbative effects due 
to poles in the expansion coefficients, and which we associate with worldline instantons. This provides 
a concrete analog of a phenomenon found recently by Drukker, Marino and Putrov in the large N 
expansion of the ABJM matrix model, in which non-perturbative effects are related to complex space- 
time instantons. In this paper we study how these very different regimes arise from an exact WKB 
analysis, and join smoothly through the magnetic region. This approach also leads to a simple proof 
of a resurgence relation found recently by Dunne and Unsal, showing that for these spectral systems 
all non-perturbative effects are subtly encoded in perturbation theory, and identifies this with the 
Picard-Fuchs equation for the quantized elliptic curve. 
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1 Introduction 

In this paper we revisit the vacuum structure oi Af = 2 supersymmetric SU (2) Yang-Mills theories 
using the resurgence formalism that unifies perturbative and non-perturbative physics^ [1-12]. The 
Af = 2 supersymmetric SU{2) Yang-Mills theory possesses a rich vacuum structure. The space of 
gauge inequivalent vacua, the moduli space, is a manifold parametrized by the scalar condensate 
u = (Tr$^). In their seminal work, Seiberg and Witten showed that this manifold is precisely the 
moduli space of genus-one Riemann surfaces (i.e. tori) and the dynamics of the low energy effective 
theory can be formulated in a geometric language in the terms of elliptic curves [13-19]. A further 
conceptual and computational breakthrough came with the introduction of the Nekrasov partition 
function [20], and the subsequent direct relation to integrable models and the Bethe ansatz [21, 22]. 
Remarkably, the prepotential in the Nekrasov-Shatashvili limit is encoded in the spectra of certain 

^Here, perturbative and non-perturbative refers to the expansion in the Nekrasov deformation parameter of the gauge 
theory, not the gauge coupling. 
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states in certain simple Schrodinger systems through the monodromy and exact WKB properties of 
differential equations [24-33]. The moduli parameter u is directly identified with the eigenvalues of 
these Schrodinger systems, and the scalar (and dual scalar) field expectation values are identified with 
actions and dual actions in an all-orders WKB analysis. 

In this paper we extend this approach to show that there are additional non-perturbative aspects of 
this relation that reflect the physics of the non-perturbatively small gaps and bands in the Schrodinger 
spectra. These spectra have three distinct physical regions, and these can be explicitly associated with 
the three physical regimes, electric, dyonic and magnetic, of the supersymmetric (SUSY) gauge theory. 
The interplay of perturbation theory and non-perturbative physics is different in each region. The 
dyonic region is characterized by divergent perturbative expansions described by resurgent trans-series 
[9, 34, 35] that systematically unify perturbative and non-perturbative physics; the electric region has 
convergent perturbative expansions but there are nevertheless non-perturbative effects associated with 
poles of the expansion coefficients. (This provides a concrete analog of a phenomenon found recently 
by Drukker, Marino and Putrov [36] in the large N expansion of the ABJM matrix model, in which 
non-perturbative effects are related to complex space-time instantons, and which were subsequently 
related to poles in the ’t Hooft expansion coefficients [37].) The magnetic region is a cross-over region 
in which non-perturbative effects are large, and in fact the spectral bands and gaps are of equal width. 
This correspondence between the Nekrasov-Shatashvili limit and monodromies of spectral problems 
also provides a simple proof of a surprising resurgence relation found by Dunne and Unsal in the 
spectra of certain quantum systems, which shows that all non-perturbative effects are subtly encoded 
in perturbation theory [9]. 

An elementary but significant observation is that the energy eigenvalue, u, for the Schrodinger 
systems should be viewed as a function of two variables, the coupling h and also the eigenvalue level 
label N: u = u{N,H). For a uniform analysis valid throughout the entire spectrum, it is natural to 
define a “’t Hooft parameter”^, 

X = Nn ( 1 . 1 ) 

and consider different limits of the two parameters, including a double-scaling limit. The usual analyses 
of divergent perturbative expansions [9, 34, 60], in the —>■ 0 limit, and their associated resurgent 
trans-series representations, are implicitly restricted to a particular non-uniform limit in which the 
eigenvalue level number N is small, N corresponding to states with energy well below the energy 
barrier. More generally, the perturbative expression for the energy eigenvalue has the form 

OO OO 

u{N, h) = u{N, A) = ^ iV2-2"i^„(A) = ^ (1.2) 

n—0 n—0 

which is analogous to the genus expansion of the free energy for a matrix model or gauge theory 
system [6]. We show that there are in fact non-perturbative corrections to this expression, of the 
form for small A, but of the form large A (see (2.30)). Moreover, in this large 

A regime the perturbative expansions are convergent, but the perturbative coefficients have poles 
that are responsible for the non-perturbative splittings. This can be compared with recent results 
concerning non-perturbative contributions to matrix models associated with the ABJM free energy 
]36-38], where in the large N limit there are extra non-perturbative terms of the form [A is 

the’t Hooft coupling], and these extra terms are related to poles of the expansion coefficients [37, 38]. 

^Note that our N is not N^. In this paper we are discussing SU{2) gauge theory. The role of Nc (or the N of the 
matrix model) is played in this context by the level number N. It is also possible to introduce yet another parameter, 
Nc, from the SU{Nc) Toda system, as in [25], but this is not done in the current paper. 
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We first review some basic facts from the SUSY gauge theory side [13-18] of the correspondence 
and set our notation, in order to make the precise identification between the gauge theory quantities 
and the corresponding spectral quantities. The vacuum expectation values of the scalar field and its 
dual partner, ao(u) and al^{u), correspond to the two independent cycles on the torus given as 

V2 r 


f V 2 

ao(u) = f /r = — / ^Ju — h? cos (j) dcj) 

J Jo 
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where A is the dynamically generated scale, and 71,72 are integration cycles discussed in detail below, 
in Section III. In this form ao{u) and Qq correspond to lowest order WKB cycles for the Mathieu 
system, hence the subscripts ”0” [14]. The information about the BPS states of the theory is also 
contained in these cycles, where the central charge and mass of a BPS state are given as 


.(m) = qeaoiu) + qmaoiu) , = V 2 \qeao{u) + qma^ {u)\. 


(1.5) 


Here qe and qm are integers that denote the electric and magnetic charge of the state. 

There are three singular points in the moduli space where one or both of the cycles acquire 
branch points. They are u —> 00 , u = and u = — A^. Around these points on the moduli space 
the low energy theory is described by weakly coupled massive Z bosons, almost massless magnetic 
monopoles and almost massless dyons, respectively. For the rest of the paper we will refer to the 
local neighborhoods around these points as the electric, magnetic and dyonic regions, respectively. 
In general, there are two separate sectors in the moduli space with different particle spectra. In one 
sector, the spectrum of the theory consists of Z bosons and an infinite tower of dyonic tower with 
charges ±(( 7 e) !)• In the other, it consists of magnetic monopoles with charge ±( 0 , 1 ), and dyons with 
±(1,±1). These two regions are separated by a closed curve, Im[a^/ao] = 0, called the curve of 
marginal stability, where the BPS states with higher charges can decay into monopoles and dyons. 
This curve is approximately an ellipse around m = 0 and contains the points u = ±A^. 

Another important object for the low energy effective theory is the prepotential, J^o(ao), which is 
a holomorphic function whose derivatives are 
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where 0ym and pym are the theta parameter and the coupling constant (at the scale oq) of the gauge 
theory. In the electric region where oq 3> A, the prepotential has the following semiclassical expansion: 
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The first two terms in this expansion are the classical and one-loop contributions, and the last term is 
the sum over non-perturbative fc-instanton corrections^ The first few terms of the instanton expansion 


comment here concerning terminology: Throughout the paper, we use the word “instanton” to describe two 
different objects: (i) the BPST instantons that appear in the 4d gauge theory; (ii) the quantum mechanical instantons 
that appear in the quantum mechanical description of the deformed gauge theory. The distinction should be apparent 
within the context of the discussion, and for the bulk of the paper it will mostly refer to (ii). 
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calculated via standard field theory methods agree with the extraction of the coefficients from the 
Seiberg- Witten solution. For example, for the Af = 2 SUSY SU{2) gauge theory, the instanton 
expansion in (1.7) is 
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However beyond the two-instanton level, the direct quantum field theory methods become computa¬ 
tionally very difficult. 

An alternative way to calculate the instanton sum is through localization [20, 21]. In this ap¬ 
proach, a two-parameter generalization of the prepotential, J^(o|ei, € 2 ), is introduced. The parameters 
characterize certain SUSY preserving space-time deformations. The prepotential of the deformed the¬ 
ory is calculable via localization technique and one obtains the Seiberg-Witten prepotential J^o(a) in 
the limit 


Jo(a) = lim T{a\ei,e 2 ). (1.9) 

ei,e2->'0 

In this paper we will focus on the Nekrasov-Shatashvili limit [21], a particular one parameter de¬ 
formation of the gauge theory with ei = li, and £2 = 0. We denote the prepotential in this limit 
as 


= lim J^(a|ei = li, £ 2 ) 

£2—>-0 

= (o, H) + (a, h) + (a, H) 

00 

= ( 1 - 10 ) 

n—0 


Continuing the example (1.8) ofN = 2 SUSY SU{2) gauge theory, the expansions for the perturbative 
and instanton contributions are 
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In the Nekrasov-Shatashvili limit, there is a direct correspondence between the gauge theory and 
integrable models [21, 22], and also with 0-1-1 dimensional (i.e. quantum mechanical) Sine-Gordon 
theory for SU{2) M = 2 SUSY gauge theory, and Lame theory for SU{2) J\f = 2* SUSY gauge theory, 
a model with a massive hypermultiplet in the adjoint representation [24-32]. In the limit where the 
mass of the hypermultiplet, m, is zero, the N = 2* theory becomes A/" = 4 theory, while at large mass 
it reduces back to the Af = 2 theory. In particular, the gauge theory moduli parameter, u, is encoded 
in the energy eigenvalue of the quantum mechanical (QM) system described by the time independent 
Schrodinger equation 


Af = 2 
Af =2* : 


2 dx'^ 
2 dx'^ 
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where V is the Weierstrass elliptic function. Some minor rescaling is required [as discussed in Section 
6.2] to have a smooth decoupling limit in which the Af = 2* theory reduces to the N = 2 theory, 
taking m? —>■ oo limit, combined with vanishing elliptic parameter k'^ —>■ 0, such that m?k‘^ is finite. 

Furthermore, at nonzero h, the periods (1.3, 1.4) generalize to the Bohr-Sommerfeld integrals of 
the QM system [25]. They can be expressed formally as all orders WKB expansions: 

OO OO 

a(u, S) = ^ a„(it) , a^{u,h) = a^{u). (1-14) 

n—0 n—0 

One recovers the Seiberg-Witten solution in the continuum limit obtained from h = 0, the leading 
order WKB approximation to (1-14). However, this is only part of the story, as these Bohr-Sommerfeld 
expressions (1-14) give only the locations of the bands or gaps. There is also important spectral 
information in the non-perturbative widths of bands and gaps. This information is connected to the 
perturbative information in quite different ways in different parts of the spectrum, corresponding to 
different semiclassical behavior for the different regions of the’t Hooft parameter A. Using the direct 
correspondence, we can associate this with different physical behavior in the corresponding SUSY 
gauge theory. 

The relation between the Nekrasov-Shatashvili limit and exact quantization conditions of quantum 
mechanical systems has also been discussed recently in [32], from the perspective of holomorphic /3- 
ensembles and the associated quantum geometry [39]. Here we follow a complementary approach, 
building our analysis on elementary WKB methods, treated uniformly across the entire spectrum. This 
relation has also been formulated in terms of Whitham dynamics [33] applied to exact quantization 
conditions for complex quantum mechanical systems. 

In this paper we build a unified WKB analysis that spans all regions of the spectrum, using 
all-orders exact WKB analysis [40-45]. In Section H we review basic properties of the Mathieu 
spectrum, and its relation to the Nekrasov partition function in the Nekrasov-Shatashvili limit. Section 
HI contains the uniform all-orders WKB analysis. In Section IV we describe a physical analog of 
the transition between different spectral regions, in terms of the transition between tunneling and 
multi-photon pair production in the Schwinger effect [46-51]. In Section V we discuss a Picard- 
Fuchs interpretation that explains the significance of resurgence, and also show that this SUSY gauge 
theory perspective in terms of the Nekrasov-Shatashvili limit yields a simple proof of the recently 
found Dunne-Unsal relation that connects the fluctuations about the perturbative vacuum with the 
fluctuations about the one-instanton sector, and with all higher multi-instanton sectors [9] . The Lame 
system is discussed in Section VI, by means of a complementary approach to the large A region, based 
on the Gelfand-Dikii expansion of the resolvent. We end with a summary and comments about future 
work. 


2 Mathieu Equation and SU{2) M = 2 SUSY Gauge Theory 


2.1 Mathieu Equation: Notation and Basic Spectral Properties 

In this Section we review relevant facts about the spectrum of the Mathieu equation [52-55] , translated 
into notation that makes explicit the relation to the Nekrasov partition function. The standard 
textbook form of the Mathieu equation is (http://dlmf.nist.gov/28): 


if” + {A —2Q cos(2t)) ip = 0 


h? d'^ip ,2 , , , 


uip 


( 2 . 1 ) 
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(We use capital letters A and Q, rather than the conventional lower-case ones, as the symbols a and 
q have special meaning in the gauge theory discussion). We thus make the identifications: 


Q = 


4A^ 

~W 



( 2 . 2 ) 


We will mostly set the scale A = 1 in what follows, re-introducing it where necessary by simple 
dimensional scaling arguments. 


m 



Figure 1. The band spectrum of the Mathieu equation, expressing the eigenvalue m as a function of the 
parameter fi, as in (2.1). This eigenvalue u will be identified with the scalar condensate moduli parameter in 
the SUSY gauge theory. The bands are shaded in grey, with the lower edges of each band shown as a solid 
(blue) line, and the top edge of each band as a dashed (red) line. At small h, the bands are exponentially 
narrow, and the band location follows the linear behavior in (2.8). At large h the gaps are exponentially narrow, 
and the gap location follows the quadratic behavior in (2.23). The top and bottom of the potential, at u = ±1, 
are shown as dotted lines. Notice the smooth transition between exponentially narrow bands (shaded) at small 
h, and exponentially narrow gaps (unshaded) at large h. This transition occurs at the top of the potential, 
where u = 1, shown as a straight line. Note that in the vicinity of the barrier top, the bands and gaps are of 
equal width, and are not exponentially narrow, as discussed Section 3.3. 

We also make explicit comparison with the work of Zinn-Justin and Jentschura [34] (see 
who used a different scaling: 

+ ^sin2(2V5x)) V' = ^ 

Thus, we identify (note: we flipped the sign of cos(a:) by a simple half-period shift) 

h=16g , u =-1 + lQg Ezjj =-1 +hEzjj (2.4) 


also [9]), 

- i)V' 
(2.3) 
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Figure 2. Another view of the band spectrum of the Mathieu equation, expressed here in terms of the 
rescaled eigenvalue E, as in (2.4), and with the scale A = 1. The rescaled eigenvalue E is to be compared 
with the normalizations in [9, 34]. The transition region between exponentially narrow bands (small h) and 
exponentially narrow gaps (large h) occurs at the top of the potential, where E = 2/h, as shown as a dotted line. 
Note that in the vicinity of the barrier top, the bands and gaps are of equal width, and are not exponentially 
small. 

Because of the direct identification of the rescaled eigenvalue u with the SUSY gauge theory scalar 
condensate, we will describe the Mathieu spectrum in terms of this eigenvalue u, with conversions to 
A 01 E made using the above re-scalings (2.2, 2.3, 2.4). 

The exact Bloch spectral condition, or Floquet analysis [52-54] , can be expressed in terms of two 
independent solutions, ipi{z;u,h) and -02(2 ;normalized at z = 0 as: 

0i(O;u, S) 'ip2iO;u,H) 

_0i(O;u,;i) u,h) 

The Bloch boundary condition, ip{z + tt) = e*® 'ipiz), can then be written in compact form in terms of 
■01 (tt; u, h) evaluated one half period away from the normalization point: 

cos (0) ='0i(7r; u, S) (2.6) 

Equation (2.6) is the “exact quantization condition”, implicitly expressing the eigenvalue u in terms 
of h [and therefore E, or A, in terms of g, or Q, using the re-scalings (2.2, 2.3, 2.4)], for each Bloch 
parameter 0. The band/gap edges correspond to 9 being an integer multiple of tt. However the exact 
quantization condition (2.6) is of limited practical use unless we have an explicit expression for the 
normalized Mathieu function 01. Concrete approximations for the eigenvalues can be obtained from 


1 0 
0 1 


(2.5) 
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the exact quantization condition by making expansions of the Mathieu functions in terms of other 
functions, such as trigonometric or Hermite functions [53-55]. Different expansions are suitable for 
different regions in the spectrum, as is familiar from elementary solid state physics. Deep inside the 
wells, we use the tight-binding approximation in terms of ‘atomic’ states bound in the wells, while far 
above the barrier we use the ‘neary-free-electron model’ [56]. 

The spectrum of the Mathieu equation consists of an infinite sequence of bands and gaps, as shown 
in Figures 1 and 2. In Figures 1 and 2, we see the transition from exponentially narrow bands low in 
the spectrum, to exponentially narrow gaps higher in the spectrum. The cross-over occurs near the 
top of the potential, where m ~ 1. In fact, this transition occurs when NH ^ ® , as discussed in Section 
3.3. In this transition region the bands and gaps are of equal width, and neither is exponentially 
narrow. In Figure I we see a transition from narrow bands, with locations approximately linear in 
N and fi., to narrow gaps, with locations approximately quadratic in N and h. Figure 2 illustrates 
the same behavior, in terms of the rescaled energy eigenvalue E = (u + l)/h, which emphasizes the 
transition from harmonic oscillator behavior at small H, to particle-on-a-circle behavior at large h. 

"electric" N f)» ^ 


"magnetic" N H ~ ^ 



Figure 3. Three different spectral regimes for the Mathieu equation: (i) far above the barrier, where u 1, 
and Nh S> 1; (ii) barrier top region, where u ~ 1, and Nh ~ 1; (hi) deep inside the wells, where u ~ —1, and 
Nh <C 1. In gauge theory language, these are the “electric”, “magnetic” and “dyonic” regions, respectively, as 
shown in the Figure. 

There are clearly three interesting spectral regimes, which we identify with three interesting phys¬ 
ical regions of the gauge theory, as shown in Figure 3: 

• dyonic: deep inside potential wells: u ~ —1, Nh <C 1, exponentially narrow bands 

• magnetic: near the barrier top: u ^ -1-1, Nh ~ 1, bands and gaps of equal width 

• electric: high above barrier top: u ^ 1, Nh ^ 1, exponentially narrow gaps 

In terms of the “ ’t Hooft coupling” (1.1), A = N h, we can consider the energy eigenvalue as a 
function of N and h, or N and A: 


u{N, h) = u{N, A) 


(2.7) 





The semiclassical limit, S —>■ 0, with A fixed, requires requires N —^ oo; but we still have three different 
regions, where A<Cl,A^l, orA^l. (In fact, we show later that the cross-over region is at 
A ~ A uniform analysis, valid for all A, therefore permits access to all regions of the spectrum. 
For example, the small A region can be interpreted as S —>■ 0, with N fixed, which therefore gives 
information about the weak-coupling expansion for low-lying energy levels with N On the other 

hand, the large A region gives information about the strong-coupling expansion for low-lying modes, 
with N fixed and 

There is an interesting inversion of the meaning of strong and weak coupling. In gauge theory 
language, the electric region is weakly coupled, while the dyonic region is strongly coupled. On the 
other hand, in quantum mechanical language in the scaling (2.3) of [34], in terms of the coupling g, 
the situation is reversed. 


2.2 Dyonic Region: Resurgent Trans-series Expansions, Deep Inside the Wells 

In this A ^ 1 regime, which is weak-coupling in the QM sense but strongly-coupled in the gauge 
theory sense, the perturbative expansions for the energy levels are divergent. The formal perturbative 
expansion reads (http://dlmf.nist.gOv/28.8.El), translated into our notation: 


u{N, h) 
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For fixed level number with <C l/Ii, the expansions for E = are factorially divergent and 
non-alternating, as series in h [9, 34]: 


E{N,g) = Y,Cn{N)n^ 

n—0 


Cu{N) 


2^^ r(n-k2Af-kl) 

T,lQn+2N+l (iV!)2 


(2.9) 


These perturbative expansions are therefore non-Borel-summable, and so are incomplete on their own, 
and should be extended to real, unambiguous trans-series expansions [9, 34]: 


fe-i 


r(iv,;i) = ^^^c„fcz(iv);i" 


n=0 k—0 l—O 


■exp(-l)' 
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In f —^ 

f^N-1/2 




( 2 . 10 ) 


The trans-series is an expansion in terms of instanton factors, each multiplied by fluctuations, and also 
by powers of logarithms due to quasi-zero modes coming from instanton/anti-instanton interactions 
[57-59]. Note that these log terms first appear at the 2-instanton level. The trans-series coefficients 
Cnki{N) are related to one another in intricate ways, such that all imaginary parts cancel, leaving a 
real and unambiguous energy [9, 34]. Mathematically speaking, the advantage of the trans-series is 
that it encodes all information about the function being computed, and is rigorously equivalent to the 
function wherever the function exists [1, 2], in contrast to an asymptotic perturbative expansion such 
as (2.8). 
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The resurgent trans-series incorporates non-perturbative contributions at all multi-instanton or¬ 
ders, the lowest of which is the exponentially small band width: 






2 24(w-ei) /2 


m 


N-l/2 


exp 



3(iV+- 




3 

+ 4 


Oih^) 


( 2 . 11 ) 


This band splitting is a single-instanton effect, and is real and unambiguous. The factor 8 in the 
exponent is the instanton action: y/2 y^cos{a;) -I- 1 dx = 8. 

However, the non-Borel-summable perturbative series in (2.8) diverges at a rate associated with 
the two-instanton sector [more precisely, the instanton/anti-instanton part thereof], and lateral Borel 
summation produces an ambiguous imaginary non-perturbative term ^ exp[—16//i]. This am¬ 
biguous imaginary non-perturbative term is in fact cancelled by an identical term coming from an 
instanton gas analysis of the instanton/anti-instanton interaction [9, 34]. This leading cancellation, at 
the two-instanton level, is just the tip of the iceberg: the cancellations between imaginary terms pro¬ 
duced by lateral Borel summation of perturbation theory, and those coming from the multi-instanton 
sectors, occur at all orders, and these cancellations are encoded in relations between the coefficients 
of the resurgent trans-series expansion. 

For example [9], for the lowest band, the large-order behavior of the perturbative coefficients is: 

c„(0) ~ n! f 1 - I • - - ^ - ... ) (2.12) 

\ 2 n 8 n{n — 1) ) 

while the fluctuations about the instanton/anti-instanton saddle are given by 

(s)‘-■ ) 

Note the precise correspondence of the coefficients in these two very different expansions. Equation 

(2.12) is associated with the fluctuations about the perturbative vacuum saddle point, while Equation 

(2.13) describes the fluctuations about the nonperturbative instanton/anti-instanton saddle point, the 
smallest action saddle-point with the same vacuum quantum numbers. This is an explicit example 
of the relations between coefficients that exist within the trans-series, and a direct manifestation of 
resurgence. In path integral language, resurgence means that the fluctuations about various different 
saddle points in the multi-instanton expansion are directly related to one another. Such resurgent 
relations persist to all orders of the non-perturbative and quasi-zero-mode expansions [9] . 

Zinn-Justin and Jentschura (ZJJ) [34] have argued the remarkable result that the entire trans-series 
can be generated from an exact quantization condition, together with just two functions: Bzjj{E, g), 
which describes the perturbative series, and Azjj{E,g), which effectively describes the fluctuations 
around the single-instanton. In fact, the resurgent trans-series structure follows naturally from a 
uniform WKB analysis, which shows it is an expression of the analytic continuation properties of the 
parabolic cylinder functions [9]. The ZJJ exact quantization condition is written as [34] (we rewrite 
the following expressions in terms ol h= IQg instead of g): 

/32\“'®^" / 32\”^^" _ 2 cos6l 

V^/ r(f-Szjj) V rd + Szjj) 
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where d is the Bloch angle, and for the Mathieu potential, the two functions Bzjj{E, h) and h) 

are given by: 




5E 
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hV/721E 525E^ 2A5E^\ 


16 J V 64 


8 


+ 


4 J 
f2^E 


- X + ^(^ + 3^') + ( 4 ; V 4 


;i\ /4487£; o 

+ 'I 6 + 


1021^;^ 


Ezn{B,h) = 1 




llE^ 


Inverting the expression for Bzjj{E,h), we obtain [9]: 


+ 3^;^ + - - 


16 


17 35E^ 25E^ 


32 


3B\ 


/ /SSB^ 205^3 405B\ 

"Viej 1^4”^ 8 '^64y“'" 


^ ^215 , 341£:2 
8 


5^4 


175^ 


(2.15) 

199E'^' 


(2.16) 


9 

32 


(2.17) 


which agrees with the perturbative expansion (2.8), with the definitions (2.4) and the identification of 
B with the level number N: 

B = N+^ (2.18) 

Thus, the function Bzjj{E,H) is equivalent to conventional Rayleigh-Schrodinger perturbation theory, 
about the perturbative vacuum. Using (2.17), the non-perturbative function Azjj{E,h) in (2.16) can 
be re-expressed as a function of B: 


Azjj{B,h) 



hV /bbB^ 205R2 135 
ley iv 4 8 '64' 


(2.19) 


The function Azjj{B,h) encodes the fluctuations about the single-instanton [9]. For example, the 
single-instanton fluctuation factor is given by 


dEzjj 

_ p 2^- 

dB 


= l_|(3i,^ + 4i,+ ^)-... 



( 2 . 20 ) 


in agreement with the fluctuation factor in (2.11). In [9], it was shown that this correspondence is 
directly connected with a simple relation between the two functions Azjj(B, H) and Ezj]{B, H): 


dEzjj 

OB 


h 

If) 


(^B + h 


dAzii \ 

dh J 


( 2 . 21 ) 


This implies that the function Azjj{B,h), and hence also Azjj{E, H), can be deduced immediately 
from knowledge of the perturbative energy Ezjj{B, h). Thus, only one of the two functions Bzjj{E, h) 
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and Az 33 {E,h) is actually needed to generate the entire trans-series. Therefore, the fluctuations 
about the single-instanton saddle, and all other non-perturbative saddles, are precisely encoded in the 
fluctuations about the perturbative vacuum. In other words, the full trans-series is encoded in the 
perturbative fluctuations around the vacuum, Ezjj{B,h). This surprising result is in fact consistent 
with the ambitious goal of resurgence, which claims that the expansion about one saddle contains, in 
principle, information about the expansions around other saddles, provided one knows how different 
saddles are connected. This connection is provided by the exact quantization condition (2.6), which is 
itself a statement of the Bloch boundary condition [9]. In Section V we discuss this further, and use 
the gauge theory perspective to give a simple proof of this surprising result (2.21). 


2.3 Electric Region: Convergent Continued-fraction Expansions, High Above the Bar¬ 
rier Top 

In the A ^ 1 regime, which is strongly-coupling in the QM sense but weakly-coupled in the gauge the¬ 
ory sense, the behavior is completely different. For small level label N, and large h, the strong-coupling 
expansions for low-lying modes, converted to our notation and with denoting the top/bottom of 
the gap, are (http://dlmf.nist.gOv/28.6.i): 
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5315625/112 

45608 


15/14 3375/18 5315625/112 


-b ... 


( 2 . 22 ) 


These expansons are in fact convergent, with a radius of convergence that increases quadratically with 
the level index N. They are generated from continued-fraction representations of the eigenvalues, 
and these continued-fraction expressions are themselves convergent. Nevertheless, despite these con¬ 
vergence properties, there are also non-perturbative effects, associated with the exponentially small 
splittings of the spectral gaps in this region of the spectrum, as are clearly seen in Figures 1 and 2. 
From a physical perspective, the expansion (2.22) governs the “fully quantum” regime where the ki¬ 
netic term dominates over the potential, and can be obtained in a straightforward fashion by treating 
the potential cos(a::) as a small perturbation to the free particle on a circle whose wave function is 
~ e^. The even/odd wave functions are identified with 6 = ttN with N being an even/odd integer. 
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Standard degenerate perturbation theory for level N then leads to (2.22) where leading term is simply 
the energy of the particle on a circle, and the higher order terms are perturbative corrections to 
it. 

Instead of taking S 3> 1 and N fixed, the high spectral region can also be probed with large A 
by taking h —>■ 0 and N —>■ oo, with ^ 1. Then for large level number N ^ 1/fi,, the continued- 
fraction expressions for the energy eigenvalues give approximate expressions for the energy of the 
gap as (http://dlmf.nist.gOv/28.6.E14): 


u{N, h) 


^ (n^ + 1 . 

9N^ + 58N‘^ + 29 \ 

'^64(fV2 _ i)5(7v 2 _ 4)(iv2 - 9) j +---I 


(2.23) 


The continued fraction relation that generates the above expansion is given in Section 3.4. Note that 
each coefficient has poles at integer values of N. In particular, the denominator of the coefficient of 
jg proportional to 


J|(fV2 _ (^2.24) 

k=l 


where [nj denotes the greatest integer less than or equal to n. 

Now we demonstrate how the energy spectrum (2.23) is identified with the multi-instanton ex¬ 
pansion of the prepotential (1.11) in the Nekrasov-Shatashvili limit of the N = 2 SU{2) SUSY gauge 
theory. First, rewrite (2.23) as: 


u 


1 [Nhy 1 / 2 y 1 5 

2 1^/ \ 64 

'a^ 1 5 1 9 1 1 2 


/ 2 y (l+ 5(w) 

\Nh) ( y/ 4^2 \ 

(Nny) (Nhy) 


+ ... (2.25) 


1 21 1 55 1 

16 a4 


+ ... (2.26) 


where we have defined the “action” 


Nh 


(2.27) 


which is half the “’t Hooft coupling” defined previously (1.1). We now compare this with the instanton 
expansion in (1.11). In order to relate the prepotential to u we use Matone’s relation [70-72], 


in OTns 


(2.28) 


The second term on the right hand side is due to the perturbative part of the prepotential. Plugging 
the expansion 


a2 a2 ^4 5^8 

A2 + 8W 25(3^ 


48ni 


vlog 


/ A4 21A8 
A2 ' ^ 256a8 


+ 


+ ... 


(2.29) 


into (2.28) one sees that the expansion for u(a, h) matches precisely with the expansion in the second 
line of (2.26). 


- 13 - 























However, it is clear that (2.23, 2.26) is not the whole story. First, this gives the same expression 
for both gap edges, whereas the expressions (2.22) clearly show that are different from one 
another. Second, the h dependence of the expansion (2.23) does not match that of (2.22). Physically, 
this is simply the statement that the perturbation theory Bohr-Sommerfeld expansion (2.23) is only 
approximate, giving the approximate location of a narrow gap high in the spectrum (see Figures 1 
and 2) for N 3> l/h, completely neglecting the non-perturbative splitting for the width of the gap. In 
fact, for a given N, the splitting of the energy levels occurs at the order l/h?^. As the level index N 
increases, the splitting drifts to higher orders in perturbation theory and becomes exponentially small 
(http://dlmf.nist.gov/28.6.E15): 




1 

4 (2^-i(//- 1)!)^ 

Nh^ / e 
27r \Nh) 



1 + 0 


A^» 1 



(2.30) 


Since the expansions (2.22) are convergent, the origin of this non-perturbative splitting is quite dif¬ 
ferent from the familiar weak-coupling analysis that associates non-perturbative terms with divergent 
perturbative series [60-62]. It is clear from the continued-fraction expansion (2.23) that the splitting 
is directly associated with the poles of the expansion coefficients at integer values of N. Physically 
this indicates the appearance of degenerate perturbation theory for both edges of the gap [56]. As 
mentioned already in the Introduction, this is a concrete analog of the results Drukker, Marino and 
Putrov [36] concerning the large N expansion of the ABJM matrix model, in which non-perturbative 
effects are related to complex space-time instantons, and which were subsequently related to poles in 
the ’t Hooft expansion coefficients [37]. We show in Section IV that in this regime the gap splitting 
(2.30) also has a natural non-perturbative interpretation in terms of semiclassical configurations. 

To make a sharper analogy with the trans-series (multi-instanton) expansion (2.10) near the 
bottom of the well, we can reorganize the large h expansion (2.22) as 


W 8 ^ f^4k ^ 8 (2^-i(V-1)!)2 ^ ;i4fc 


(2.31) 


where the coefficients ak{N) are those given in the continued fraction expansion (2.23). However 
this identification is valid only up to order k = {N — l)/2, due to the poles in the coefficients of the 
continued fraction expansion. The gap splitting, encoded in the PkiN) terms, arises beyond this order, 
and has leading behavior given in (2.30). So, at finite N the gap splitting terms are missed by the 
Bohr-Sommerfeld expansion. 

In the next Section we show how we can actually obtain these non-perturbative gap splittings 
from an exact all-orders WKB analysis. In Section 4 we will further identify this effect with worldline 
instantons, and make a direct physical analogy with multi-photon ionization in monochromatic time- 
dependent laser pulses. 


3 All-orders WKB Analysis of the Mathieu Equation: Actions and Dual 
Actions 


The all-orders WKB expression for the location of bands and gaps can be expressed as [40-45] : 

2ttN h (gap) 


P = 


2t: {N ■ 


(band) 


(3.1) 
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where P is the local momentum and the contours go around the appropriate turning points, as shown 
in Figure 4. Geometrically, for the Mathieu problem there are two independent cycles and they 
correspond to the generators of the two cycles of the torus. More explicitly, with u denoting the 
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Cm ■ ■ ^ 
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-2tt 

0 * : 
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Cg 



-iTT 
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* t 



Figure 4. The integration contours for the all-orders WKB expressions. The first figure shows the paths for 
the bands, when the energy is well below the barrier top. The contour Cb is used to determine the location 
of the band, while Cb is used to determine the width of the band. The second figure shows the paths for the 
gaps, when the energy is well above the barrier top. The contour Cg is used to determine the location of the 
gap, while Cg is used to determine the width of the gap. Note that in this case the relevant turning points are 
in the complex plane. 


energy eigenvalue, we have [40, 41] 


1 

27r 



27r 


Vu 


Vdx — 

26 


/ jvr , 


Nh (gap) 

[N + (band) 


49(G')^ 

(u-l/)ii/2 


IQV'V'” \ , 

{U - 1 /) 7/2 ) 

(3.2) 


With proper analytic continuation, this quantization condition permits smooth transitions and duali¬ 
ties between the various spectral regions, connecting weak and strong coupling, and also the bottom 
and top of the wells. The distinction between the various regions is encoded in the location of the 
turning points in the complex plane, and the associated Stokes lines [44]. For energies inside the 
wells there are real turning points. As the energy approaches the barrier top, the turning points come 
together and coalesce, and move apart again along the imaginary axis for energy above the barrier 
top. See Figures 4 and 5 . 

For example, the leading behavior in (2.8) and (2.23) can be found immediately as follows. For 
the location of the center of a high gap: 

^~ ° (i) 

For the location of the center of a deep band: 

(^ + 2 ) ^ ~ ~ ]lI Zy2 ~ (1 + m) + O ((1 -b u)^) 

(3.4) 


u{N,h) -^ + -.- (3.3) 


h) —1 h ^ N — ) -b... 


The higher terms are discussed below, in the next subsections. 

Perhaps less well-known is that there are corresponding expressions for the width of a band deep 
in the spectrum and of a gap high in the spectrum [63, 64, 66-68]. At leading order, the width of a 
band or a gap is expressed as the product of a density-of-states factor and an exponential: 


Autv 


2 du 


-Im 

h 



(3.5) 
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Figure 5. To determine the location of the bands or gaps, we integrate the WKB integrands along the paths, 
Tband and Lgap, shown here. In the deep band case, the width of the band is determined by integrating through 
the barrier. As the energy goes above the barrier top, the turning points coalesce and move into the complex 
plane: in the high gap case, the width of the gap is determined by integrating between these complex turning 
points. 


where now the contours are around the dual contours Cb or Cg shown in Figure 4. Leading order 
results for bands and gaps in the barrier region are discussed in [65-68], and in more detail below in 
Section 3.3. 

As discussed in the introduction, the WKB periods associated with the Mathieu and Lame equa¬ 
tions correspond to the vacuum expectation value of the scalar field $ and its dual partner, and the 
energy eigenvalue u correspond to the gauge invariant observable (Tr$^). We set our normalization 
such that classically $ci(a) = |(t^, where cr^ is the Pauli matrix, and therefore a = -h ..., 
fixing the normalization to 


a(u) 


aoiu) 


LI P = 





(W) 


/^2 


2 ^cos~^(ii) 


dx — . 


/ cos J^2 pcos 

y/u — V dx — 

-cos~^(u) ^ J — cos~^ 


(W) 


/^2 


(u) 


(u - F)5/2 


dx — 


(3.6) 

(3.7) 


3.1 Dyonic Region: Resurgence from All-Orders WKB 

For energies below the barrier top, where —1 < u < 1, the leading WKB expression for the band 
location is given by the usual Bohr-Sommerfeld expression that involves only the real part of the 
action, while the dual action is associated with under-the-barrier, and so is pure imaginary^. The real 
action in this region can be expressed as a linear combination of the two independent actions defined 
above: 


7le[a{u)] = a{u) + a^{u) 


(3.8) 


The exact quantization that identifies u with the center of the band is implemented by requiring 


Re [a(u, h)] = a{u) + a^{u) 



(3.9) 


'^Here, and for the rest of the paper, the notion of “rear’/“pure imaginary” action refers to the particular choice of 
cycles such that when —1 < u < oo, and h £ R"*", the associated action is real/pure imaginary. In general, u and h can 
take complex values, and the analysis follows by analytic continuation. 
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The leading order terms of the actions are expressed in terms of elliptic integrals 


Tle[ao{u)] = 


V2 


dy 


y-u 


y^ — 1 IT 


E 


M + 1 




It + 1 


Df \ 

oo (w) =- 


dy 




2 / - 1 


1 — u 


- 2 (^+ 1 )^ 


1 — u 


(3.10) 

(3.11) 


These expressions are the original Seiberg-Witten solution, which in the quantum mechanical language 
is identified with the leading order WKB expansion. Note that they satisfy the Picard-Fuchs relation: 


ao(u) 


dttQ {u) D. .dao{u) 2i 
— Gq (u)- — 


(3.12) 


du du TT 

which follows from the Legendre relation [here K' denotes K(1 — fc^), using the conventions of [52]]: 

EK'+ E'K - KK'= I. (3.13) 


The Picard-Fuchs equation (3.12) is invariant under S'L(2,Z) transformations applied to the pair 
(ao,a^) which corresponds to changing the basis for the two cycles [69-72]. For example, we can 
replace Gq by i?e[ao] in (3.12) and the equation equation will still hold. 

Remarkably, the higher-order WKB actions can be obtained by acting on these leading WKB ac¬ 
tions with differential operators with respect to the energy u [25, 27]. This follows from the fact that 
for V = cos(a::), the numerators in (3.7), which are given by the derivatives of V, can be re-expressed 
as polynomials of V. Therefore by differentiating \/u—V with respect to u taking appropriate combi¬ 
nations one can generate the integrands in (3.7) up to total derivatives which vanish after integrating 
over turning points. For example, at the next two orders: 


, , ^ d^ d\ 


( 3 . 14 ) 

( 3 . 15 ) 


and the same relations hold for g^{u) as well. The general form of the differential operators that 
relates a„ to Gq is 


n 

an{u) = 


(n) fc d"+'^ao(a) 
k du'^+^ 


(3.16) 


We have verified (3.16) to the order h}°. The coefficients up to this order are given in Appendix 
B. The first two next-to-leading order actions calculated from (3.15) are 

= 48x((-..=) (<' “ (^) + (^)) 
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Tle[a2{u)\ = - 


02 (u) = 


4608071 (1 - v?y 


4608071 (1 - v?y 


(1 - u){4u^ + 930^ - 60u + 75)K 


1 + w 


+2 - 153u^ - 75) E 


1 + M 


(1 + u)(-4it3 + 930^ - 60u^ + 75)IK 


1 — u 


+2 (4u^ - 153u^ - 75) E 


1 — u 


(3.19) 


(3.20) 


The actions can be expanded about the bottom of the wells, u ~ —1, as: 


Re [ao(u)] 
Re [ai(u)] 
Re [ 02 ( 0 )] 


M + 1 {u + iy 3{u + iy 25{u + iy 245(0 + 1)^ 

2 32 512 16384 524288 ' 

1 5(m + 1) 35(u + 1)2 525(0 + 1)3 8085(0 + 1)^ 

m 2048 32768 1048576 33554432 

17 721(u + l) 10941(0 + 1)2 141757(0+ 1)3 

262144 8388608 134217728 2147483648 ' 


(3.21) 

(3.22) 

(3.23) 


Recalling (2.4) that 1 + o = HE, we can re-write these expressions for 77,e[a„(o)] as expansions in 
powers of h and E: 
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. (3.24) 

(3.25) 

(3.26) 


Comparing with (2.15), we recognize these expansions as the highest powers of E, the next-to-highest 
powers of E, and next-to-next-to-highest powers of E, respectively, for each power of h in the function 
Bzjj{E, H). Thus we have the following identification with the results of Zinn-Justin and Jentschura 
[34]: 

2 ^ 

Bzjj{E, fi) = K(-l + f^E)] (3.27) 

n—0 

= ^Re[a{-l + hE,h)] (3.28) 

Thus the “exact Bohr-Sommerfeld condition” (3.9) for the location of the energy bands deep inside 
the wells expresses the perturbative expansion for the location of the band. Inverting this Bohr- 
Sommerfeld condition, to express u = u{N, h) as a function of N and fi., we arrive at the perturbative 
expansion (2.8), or in the notation of ZJJ, the expression (2.17) for the energy E{B, H) as an expansion 
in h and B = N + ^. Recall that this expansion (2.17) is non-Borel-summable. 
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Similarly, the dual actions can be expanded about the bottom of the wells, m ~ — 1, as: 
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Notice that for n > 1, the {u) diverge as u —)■ —1. 

The dual action a^{u,h) is related to Zinn-Justin’s Azjj{E,g) function as follows: 

Azjj{E,h) = ^a^{-l + hE)-2\nr(^+B{E,h)) +H2 tt) - 2B{E,H) \n 


(3.30) 


(3.31) 


(3.32) 


= ~ Im{a^{-l + hE)) - 21nr (^^+B{E,h)^ + ln(27r) - 2B{E,h) In 


(3.33) 


The comparison for A(E,h) is non-trivial, as it requires using the large B asymptotics of the InT 
function. The subtraction of these terms corresponds to matching the perturbative solution near the 
bottom of the well to the solution near the top of the barrier, coming from the exact solution for an 
inverted harmonic well, which is how a^{u) looks near the top of the barrier [66, 67]. 


3.2 Electric Region: Convergent Expansions and the Nekrasov Instanton Expansion 

In the electric region (i.e. u —>■ oo) the real period is identified with a(u) and the pure imaginary 
period is still a^{u). The exact quantization that identifies u with the center of the gap is 


a(u) = 


h 
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N 


The actions can be expanded for m >> 1 as: 


,,/—/! 15 105 \ 
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(3.34) 


(3.35) 

(3.36) 

(3.37) 
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Combining these expansions we find 
2 

- {ao{u) + h^ai{u) + fi^a2{u) + ...) 
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16m2 1024w4 16384m6 ''' 

^ 35 1155 75075 

8(2lt)5/2 

32w2 1024 m4 65536m6 

273 5005 2297295 
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32(2m)7/2 

' 64m2 512m4 131072^6 


+ ... 


(3.38) 


Identifying the left-hand-side with N^ and inverting, we obtain the expansion (2.26) of the gap energy 
u{N, h). Thus, the all-orders-WKB action a(u, h) determines the (convergent) expansion of the location 
of the gap high up in the spectrum. 

The dual actions can be expanded for m 3> 1 as: 
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—wf* (u) 
-wf (u) 




1 


-2 -I- log(8M) 


l-log(8u) 47- 301og(8u) 


16it^ 


2048^4 


-f .. 


2471-^/^ 

1 


13 - 61og(8M) 883 - 4201og(8u) 


452®7r(2M)^/^ 


-1-f 


16u2 1024^4 

567 - 180 log(8M) 127461 - 49140 log(8M) 
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1024^4 


(3.39) 

(3.40) 

(3.41) 


These dual actions determine the exponentially narrow width of the gap, high up in the spectrum, as 
discussed in the next subsection. 


3.3 Magnetic Region: Duality and Analytic Continuation Across the Barrier 


Across the magnetic region, there is a transition from the divergent perturbative behavior character¬ 
istic of the dyonic region, and the convergent perturbative expansions characteristic of the electric 
region. As is clear from the plots in Figures 1 and 2, the transition is smooth, but connecting the 
regions requires careful interpretation of the various expansions. Of particular interest are the different 
mechanisms by which non-perturbative terms arise in the different physical regions. For example, the 
general expression for the exponentially narrow width of a band deep in the dyonic region, and of a 
gap high in the electric region is [66-68]: 
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(3.42) 


In the dyonic region, 
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(3.43) 

(3.44) 


Therefore, from (3.42) we obtain the band width estimate (using Stirling’s formula in the last step): 
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(3.45) 
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in agreement with (2.11). On the other hand, in the electric region, 

1 2 t2 

^ 2^0 ■ ■ 

7 rlm[a^] ^ (ln(8u) — 2) + ... 

Therefore, from (3.42) we obtain the gap width estimate: 




h^N / e \ 
27 t V HN ) 


2N 


(3.46) 

(3.47) 

(3.48) 


in agreement with (2.30). Thus, the formula (3.42) has the correct form in both extreme limits, in one 
case referring to the width of a band, and in the other to the width of a gap. 

The magnetic region near u ~ 1 is more subtle. In this regime, 
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u — 1 
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(3.49) 

(3.50) 


The latter relation tells us that the exponential behavior becomes of order unity. The first relation 
gives us the leading scaling between N and H: 
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(3.51) 


It is clear that the barrier top is in the vicinity of fV ~ l/fi., but the above fixes the non-trivial 
coefficient to be S/tt. It is instructive to evaluate the energy eigenvalue u{N,h) with this scaling, in 
both the dyonic and electric regions, using (2.8) and (2.23) respectively: 
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(3.53) 


It is remarkable that these two very different expansions coincide at it ^ 1. 

In fact, we can refine further the estimate in (3.51). The edges of the bands/gaps when it = 1 are 
given by [68] 


4 ttS 


(3.54) 


as shown in Figure 6. We also see that in the immediate vicinity of these points, the dependence of 
It on Q{= ^) is approximately linear, with a slope depending inversely quadratically on TV: 
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(3.55) 


This implies that in the vicinity of the barrier top, where it « 1, the bands and gaps are of equal width, 
and are not exponentially narrow. This can be seen clearly in Figure 6, where from a given intersection 
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Figure 6. Plots of the bands (shaded) and gaps (unshaded) in the magnetic region, in the vicinity of the 
barrier top where u = 1. The vertical lines mark the values Q = ^ (^N ± where N is the band label, and 
which coincide very accurately with the points at which the band/gap edges intersect the line u = 1. 


point at M = 1, we observe equal widths of the band and gap above and below that intersection point. 
In fact, 


~ ~ 0{h) 


(3.56) 


Another way to understand the smooth transition across the barrier top is to note that in both 
the dyonic and electric region, we must avoid poles by analytically continuing h (or equivalently g) 
off the positive real axis, which effectively gives a small imaginary part to u. This then avoids the 
divergences of an(u) in the region u ^ 1, in a way that connects smoothly. Since all the actions are 
expressed in terms of the elliptic functions E and K, this relies on the analytic continuation properties 
of these functions. More precisely, the analytical continuation properties connect different points in 
the moduli space as follows 


/ ao(-M) 0\/ ao{u) \ 

yaE{-u)) - * U ±v w(«)y 


(3.57) 


where ±1 denotes sign{X7n[u\). This is the manifestation of the residual Z 2 symmetry of the broken 
U{l)]i symmetry of the gauge theory [71]. For the subleading terms from (3.16) we see that 



-(- 1 )" 


/±1 0 \ /a„(u)\ 

Iv 1 ±lj {a^iu)J ■ 


(3.58) 


In particular, for —1 < m < 1, this relation can be used to relate the magnetic regime {u ^ 1) to the 
dyonic regime {u ^ —1). The appearance of sign{Xm[u\) reflects the necessity of analytic continuation 
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mentioned above, and is ultimately related with the fact that h G is a Stokes line for the small h 
expansion. 

Also, the region with 0 < u < 1, can be related to the 1 < m < oo region by invoking the inversion 
relations for the elliptic functions E and K. These elliptic functions have non-trivial connection 
formulas, with ambiguous non-perturbative imaginary contributions when evaluated on the real line. 
See, for example, http://dlmf.nist.gOv/19.7.E3 
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(3.60) 


where ± is determined by the sign of Im{k'^).^ With the above connection formulas, the transition 
across the barrier is smooth. 


3.4 Connecting Strong and Weak Coupling Regimes 

In this section we explain the connection between the large h and small h expansions. These two 
expansions are very different in nature: the former is convergent and the latter is divergent. Yet 
the exact WKB methods described in the previous section can be used to relate one to the other. 
We present an explicit method for generating the small h expansion for the low lying modes (i.e. 
A ^ 1), which governs the dyonic region, using the knowledge acquired from the “continued fraction 
expansion” for A 3> 1, which governs the electric region. 

In the electric region, the energy eigenvalue u for high lying states can be obtained from the 
continued fraction expression 
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(3.61) 


as an expansion in = 16/S^, by making a series ansatz for u and equating both sides of (3.61). 
This expansion is valid for high lying levels with N ^ 1/h. The hrst few terms of it are given in 
(2.23). Next, we substitute a for N by using the relation a = Nh/2, which is the exact quantization 
condition (3.34) in the electric region. In the electric region a = A/2 ^ 1 and we reorganize u as a 
series in Xja^ and h, as written in (2.26). Having obtained this double series, we invert it and obtain 
a as a double series in u and H, i.e. a{u,h) = fe which is valid for u S> 1. So 

far, this is an alternative way of obtaining the large A expansion (3.38), directly from the continued 
fraction expression (3.61), complementary to the WKB derivation decribed earlier in this section. 

However, from the general form of the WKB period a(u,H) given in (3.16), we know that each 
coefficient of is related to the lowest order coefficient [with n = 0] via a differential operator, with 
expansion coefficients Therefore, for any given order by applying the differential operator 

in (3.16) to the coefficient of hP of (3.38), and comparing with the coefficient of we can solve for 
K^"^s. Once we know these coefficients, we can construct both a^{u) and a{u) from a^iu) and 
Oo(u), for any ul 

To obtain the small h, small N (dyonic) expansion for u, we then expand a and around u ~ — 1. 
In the dyonic region, the exact quantization condition is 'lZe[a] = a + = h/2{N + 1/2). Then we 


^Note that these analytic continuation properties are incorrectly stated, without the ± signs, in many books and 
tables [73]. 
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invert back to obtain u(TZe[a\,H) which is the perturbative, small h expansion. The non-perturbative, 
instanton induced part of the trans-series can be obtained from as described in Section 3.1. 

Before closing this section, we also would like to point out that the continued fraction expansion of 
u{a, h) given in (2.23) can alternatively be obtained from the conformal block expansion obtained from 
the AGT correspondence [23, 24] in the £2 —> 0 limit. The instanton part of the Nekrasov partition 
function® has a rather simple group theoretical expansion given as [25, 74] 

/a ' \ ^ / \2 \ 

Z“t(a;ei,e2)=exp .F(a;ei,e2)) =5] -) QX'([l”].r]) (3-62) 

V ere2 / „ 

^ n—0 ^ 


where QA{y,Y') = (A|LyL_v/|A) is the Shapovalov matrix associated with a conformal primary 
|A) with Young tableaux Y and Y' satisfying |y| = |Y'|, and [1"] is a shorthand notation for Y = 
{1,1,..., 1} with |Y| = n. The first couple of relevant entries of the inverse of Qa are 


Qa([1]J1]) = ^ , QKX[11],[11]) = 


8 A- 


4A(16A2 -h 2cA - lOA c) 


(3.63) 


The higher order terms can be computed in a straightforward way. The AGT correspondence maps 
the conformal dimension, A, and central charge, c, to the gauge theory variables: 


A = — , c = 1 - . 

£162 V 4 / £162 

In the £2 —> 0 limit, the prepotential stays finite: 

J-{ri‘ (a; £i) = 1 ™ £2 log £ 1 , £ 2 )) 
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(3.65) 


Using Matone’s relation (2.28) and switching from a to N using the exact quantization we get 
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(3.66) 


which exactly coincides with the expansion (2.23) we obtained from the continued fraction (3.61). The 
leading term, N'^, comes from the perturbative part of the Nekrasov expansion. 


4 Worldline Instantons and Multi-photon Vacuum Pair Production 

As discussed in Section 3, despite the fact that the perturbative expansions in the electric region are 
convergent, there are still non-perturbative instanton effects, leading to narrow gap splittings high in 
the spectrum. These splittings (3.48) have a very different form from those (3.45) low in the spectrum. 
This is analogous to the results of Drukker, Marino and Putrov [36] for the large N expansion of the 
ABJM matrix model, in which non-perturbative effects are related to complex space-time instantons. 
Here we give an analogous physical identification with worldline instantons [75], for this effect in the 
Mathieu spectrum. Furthermore, we present a simple physical interpretation of the transition between 
different spectral regions, in terms of the transition from tunneling pair production to multi-photon 

®Note that our overall normalization of tF differs from some of the literature on the subject. 
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pair production. This analogy gives a physical example in which a non-perturbative quantity, the pair 
production probability, turns smoothly into a multi-photon expression of a very different form. 

Recall the semiclassical result of Brezin and Itzykson [50], and Popov [51], (generalizing Keldysh’s 
work in atomic ionization [46]), that for a monochromatic time dependent electric held, £(t) = 
£ cos(a;t), the usual Schwinger formula [47-49] for the vacuum pair production probability in a back¬ 
ground electric held becomes 


probability 


exp 


m? TT 


5 ( 7 ) 


(4.1) 


Here the dimensionless “Keldysh adiabaticity parameter” is defined as 

rriLu 



(4.2) 


For this monochromatic time-dependent held £{t) = £ cos(a;t), one hnds for the function ( 7 ( 7 ) in (4.1) 
the expression [49-51, 75]: 
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In the static limit, 7 —>■ 0, we recover the familiar Heisenberg-Schwinger result, 

r 2 

Tfl^ 

probability ^ exp- — 


(4.3) 

(4.4) 


(4.5) 


which is obviously non-perturbative in the strength of the applied held, and has a well-known inter¬ 
pretation in terms of tunneling from the Dirac sea through the barrier created by the constant electric 
held. In the opposite limit of a high frequency held, where 7 ^ 1 , the logarithmic beahvior of 3 ( 7 ), 
shown in (4.4), implies that the semiclassical result (4.1) in fact produces a perturbative result: 


probability 


exp 
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(4.6) 


The physical interpretation of (4.6) is that ‘Imjijj is the multi-photon number, the number of photons 
of energy w required to excite the virtual pair over the “binding energy” 2 mc^, by a multi-photon 
process, rather than by tunneling [46, 50, 51]. It’s a probability, so the power is twice this multi¬ 
photon number. The hnal answer is perturbative, as it is a power of the applied electric held, but it is 
a very high power in the relevant semiclassical limit where w <§1 m. Thus we see that the semiclassical 
expression (4.1) has two very different limits, and we show below that this is analogous to the transition 
from exponentially narrow bands (3.45) low in the Mathieu spectrum, to power-law narrow gaps (3.48) 
high in the Mathieu spectrum. 

First, recall that the above pair production results have a natural interpretation in terms of 
worldline instantons [75], which are hnite action periodic solutions to the classical Euclidean equations 
of motion for the world line of a charged particle in the background electric field: 




(4.7) 
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For such worldline instanton solutions, a semiclassical approximation to the Feynman’s world line path 
integral representation of the QED effective action gives the leading expression for the pair production 
probability 

probability -- exp [-5[a;i„stanton]] , S'[x] = J dr (4.8) 

For the cosine electric field, £{t) = S cos{ujt), the classical trajectories are known, and the associated 
worldline instanton action indeed equals the exponent in (4.1, 4.3) [75]. 

The relation to the Mathieu equation is the following. The pair production probability can alter¬ 
natively be computed via a Bogoliubov transformation as a quantum mechanical reflection coefficient 
in the Klein-Gordon equation (we take zero electron/positron momentum, to get the leading effect, 
and treat scalar QED instead of spinor QED, since this also gives the same exponential behavior) 
[50, 51]: 


— (j) — [ -I- ( — sin(w t) 
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cos(2x)(^ =(409) 


written now in Mathieu form. So we identify the Mathieu equation parameters: 


A = 






Therefore, we further identify [using the scalings in (2.2) and (2.4)] 


2 _4uj^ 




Further, note that the leading dual action in (3.11) can be re-expressed as 


Therefore, 


This implies that^ 
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(4.15) 


To see how this works in the various limits, recall that in the semiclassical approach to vacuum pair 
production, the electron mass sets the dominant scale, so we require 8 <C and to m. But this 

^Note that the pair production probability is related to a bounce, which is a closed-path instanton/anti-instanton 
configuration, referred to as a “worldline instanton”; while the Mathieu band or gap splitting is a single-instanton effect. 
For time-dependent electric fields associated with more realistic laser pulses, there are important quantum interference 
effects and these are captured by complex instanton trajectories [76]. 
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still permits arbitrary values of the adiabaticity parameter 7 = mw/£. In the static limit, 7 1, this 

implies w ^ 1 , and 


dTT ^ r n/ XT 27r , _ 47r7^ ttto^ ,. , 

— (4.16) 

which leads to the standard Schwinger formula for pair production in a static electric field. 

On the other hand, in the multi-photon limit, 7 ^ 1 , we have 

- » 4 (“-U) 

so we can consistently consider the hierarchy of scales 

£ moj ^ w? (4.18) 


Then defining N = ^, we see that it corresponds to the gap label index N in the Mathieu spectral 
problem. Indeed, in this limit 


rn £ A 
N = — ^ ^ = - 


(jj 




(4.19) 


in Mathieu language. Thus, from the Mathieu equation gap splitting formula (2.30) for large N [using 
Stirling’s formula], the identification (4.11) leads to 


(gap width) ^ 



(4.20) 


which, up to a multiplicative pre-factor, is the Brezin/Itzykson/Popov multi-photon expression (4.6). 


5 


A Simple Proof of the Dunne-Unsal Relation and its Geometric Inter¬ 
pretation 


In this Section we present a proof of the Dunne-Unsal relation (2.21) from the gauge theory point 
of view. In terms of quantum mechanics, this relation remarkably links the perturbative fluctuations 
around the vacuum to the non-perturbative fluctuations around instantons [9] , as mentioned in Section 
2.3. In addition to the SUSY inspired proof, we also show that the Dunne-Unsal relation can be 
identified as the generalization of the Picard-Fuchs equation for the quantized elliptic curve with 
nonzero h. Physically this identification is an explicit example of the connection between the resurgent 
trans-series expansion and the geometry of compact Riemann surfaces. 

Our starting point of the proof is the generalization of Matone’s relation for h ^ 0, 


u{a, H) 


ZTT ^9J^Ars(a, S) 

dA ~ Y, ■ 


(5.1) 


where as mentioned in Section 2.3, the instanton part of the original relation [70] is unchanged and 
the shift in u is due to the perturbative contribution to [33]. The next step is to observe 

that the the prepotential can be expressed as follows 


J^Ns{a, K) := 
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(5.2) 


- 27 - 






where the we use the symbol" to denote dimensionless quantities.® This follows from rescaling the 
Schrodinger equation (1.13) and the corresponding periods (3.7). After the rescaling, it follows that 






dA da dh 

OT- t dTNs{0',h) dTNs{a,li) 

= 2T„,(a.H)-a - - -fi- - - 

Then, by differentiating (5.1) with respect to a and using (5.4) we obtain 


2idu(a,h) Tj da^{a,h) da^(a,h) 

TT da da da 


(5.3) 

(5.4) 


(5.5) 


where we used the fact that dT/da = . The hnal step in order to arrive at the form which relates the 

instanton expansion to perturbative expansion is to switch to the variables Azjj{B, H) (which encodes 
the fluctuations around instantons) and Ezjj{B, h) (which encodes the perturbation expansion) where 
Ezjj = {u— l)//i, B = 2a/h, and Azjj is deflned as in (3.33). In these variables, (5.5) becomes 
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dB 
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16 
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5Azjj \ 
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It is also illuminating to switch the independent variables from (a, S) to {u,h). With this change 
of variables, the terms in (5.5) transform as follows: 

du_/daY^ da^ _ 9 a^ dada’^^ 

l)a~\du) ’ Ihj^ ’ ^lm~\du) \lmdu~ dhlhj^) 


where the independent variables on the left-hand side are {a,h) and on the right-hand side are (u,h). 
Then, the Dunne -Unsal relation takes the form of Picard-Fuchs equation, extended to nonzero h: 



daP 

du 
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(5.8) 


Note that (5.8) is invariant under the SL{2, Z) transformations that act on the pair (a, a^) as expected. 
It is further useful to express (5.8) in terms of the expansions (1.14): 
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(5.9) 

(5.10) 


The second line of this equation shows that all the higher order terms contribute as zero to the Picard- 
Fuchs equation, and the constant 2i/Tr on the right hand side does not get any corrections at nonzero 
h. In other words, by using the Riemann bilinear identities on (5.10) and (5.10) we deduce that only 
the zeroth order periods oq, a// contribute to the “total flux” on the torus, whereas the higher order 
terms do not. However, the cancellation of the flux from higher order terms is a result of rather 
intricate cancellations that involve contributions from different orders as seen in (5.10). 

The identification of the Dunne-Unsal relation with the Picard-Fuchs equation provides an explicit 
geometric interpretation of the connection between perturbative series and instanton expansion. The 

®Note that h has mass dimension 1 in the gauge theory/QM correspondence. This can easily be seen from the 
Schrodinger equation (1.13). 
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perturbative corrections to the energy eigenvalue are characterized by the quantization of the real 
period, while the exponentially suppressed instanton corrections are characterized by the dual pure 
imaginary period. The Picard-Fuchs relation connects these two independent periods. Furthermore, 
even though the perturbative and non-perturbative expansions take very different forms in different 
regions of the energy spectrum, this connection holds throughout the spectrum. 

It is interesting to note that similar connection formulas between perturbative and non-perturbative 
physics exist for other QM potentials, such as the double-well and SUSY double-well [9]. In fact, in 
general a relation like (2.21) exists whenever the spectrum is such that the associated Riemann surface 
is described by just the two dual actions, a and . For more general cases, where the Riemann surface 
is of higher genus, the relation will generalize to include multiple pairs of actions whose number is 
equal to the genus [77]. 


6 Lame Equation and SU{2) Af = 2* SUSY Gauge Theory 

As discussed in Section 3 for the Mathieu equation, a formal expansion of the action may be obtained 
by an all-orders WKB analysis. This can then be expanded in the high or low energy region. Instead 
of repeating these steps for the Lame system, we summarize the novel features that arise in the dyonic 
region in Section 6.1, and then use a completely different technique to study the electric region, using 
the relation between WKB and the KdV hierarchy: see Section 6.2. 


6.1 Resurgent Analysis in the Dyonic Region 

The dyonic gauge region of the N = 2* theory corresponds to the A ^ 1 spectral region of the 
Lame system. Previous analyses [27] have not taken into account the fact that in this region the 
spectral expansions are divergent and non-Borel-summable, as for the Mathieu system, and so should 
be described by a resurgent trans-series. In fact, the associated resurgent structure is even richer 
than for the Mathieu system, due to the existence of both real and complex instantons [35] . This is 
ultimately due to the fact that the Lame potential is doubly-periodic in the complex plane, and even 
though the quantum mechanical path integral is a sum over real path configurations, the existence of 
complex saddle paths has a direct influence on the divergent structure of perturbation theory [35] . 
The Jacobian form of the Lame equation is conventionally written http://dlmf.nist.gov/29. 

2 .i: 

+ {H — v{v + l)/c^sn^ (z, = 0 (6.1) 

The Jacobi elliptic function sn^ is related to the Weierstrass 7^-function as®: 

r {x + iK';K,iK') = - \ {k^ + 1) + fc® sn®(a:; /c®) (6.2) 

o 

where K' = IK(1 — fc®). The Lame equation (6.1) reduces to the Mathieu system in a special scaling 
limit: since sn®(z; 0 ) = sin®(z), we must combine the fc® —>■ 0 limit with the v ^ oo limit, in such a way 
that the combination k® = viy -\- 1) fc® remains finite. This then leads to the natural identifications: 
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\Jv(v + l)k‘^ 


u —1 + — H 


(6.3) 


^We use the standard period conventions of [52]; other papers sometimes use differing conventions [27, 30, 88, 89]. 
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With these identifications, we can express the standard perturbative expansions of the Lame eigenval¬ 
ues (http://dlmf.nist.gOv/29.7.El) as: 


u{N, h) = —i-\-h 
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28 







(6.4) 


Notice that we recover the Mathieu expression (2.8) in the limit —>■ 0. 

For a given N <^l/h and these expansions are divergent, with factorially growing coefficients. 
But the large order behavior can be alternating or non-alternating, depending on This more 
intricate structure is due to the existence of both real and complex instantons, and is discussed 
in detail in [35].^° In particular, it means that the leading large-order growth of the perturbative 
coefficients is not solely governed by the real instanton/anti-instanton action, but also by the complex 
(‘ghost’) instanton/anti-instanton action. 

Associated with this resurgent divergent structure of the perturbative expansions is the existence 
of non-perturbative band splittings, for any . Approximate expressions for these band splittings are 
given in (http://dlmf.nist.gOv/29.7.E5); for details see [78, 79]. Translating these results to our 
notation (6.3), the splitting of the A^*'' band is 

~ (6.5) 

2h [2 r 32 

mV ^ \H1 - fc2 )) -pfc j j 







This band-splitting is a one-instanton effect, as the instanton action for the Lame potential is [80, 81] 





( 6 . 6 ) 


Note that this reduces to the Mathieu expression (2.11) as —>■ 0, recalling that (1 — k)^/^ —>■ 1/e in 

this limit. 


6.2 WKB, Gelfand-Dikii Expansion, and KdV, in the Electric Region 

As discussed in Section 3 for the Mathieu equation, the electric region of the Lame spectrum can be 
studied by using a large u expansion of the all-orders WKB actions a{u,h) and a^{u,h). Here we 
describe a different, complementary, technique to analyze the electric region, using the work of Gelfand 
and Dikii concerning the high-energy asymptotic expansion of the resolvent [82, 83]. The coefficients 
of this expansion are expressed directly in terms of the KdV integrals, evaluated on the QM potential. 
This observation is particularly interesting for the Lame equation because the expansion coefficients are 
simple polynomials, with interesting combinatorial properties, in the scaling parameter that multiplies 
the elliptic potential [84]. This is described below. 

For the Schrodinger equation with a periodic potential V(x) of period L, 


dx'^ 


tp + V(x) iIj = E tp 


(6.7) 


^^The paper [35] uses another form of the elliptic potential, sd^{z,k‘^), which has a manifest symmetry under k? 
1 — k^. This potential is related to the sn^( 2 :, k"^) by a simple Landen transformation. 
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the Gelfand-Dikii expansion [82, 83] relates the level number N to the high energy asymptotics as 


ttN j — 1 —A 

VE-- 1 : 


ij+i 
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( 6 . 8 ) 


where Ij+i\y\ are functionals of the potential V{x) given by the KdV conserved quantities. These can 
be generated by simple recursion relations, and the first few are: 



f Vdx 

0 

pL 

(6.9) 
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(6.10) 

h[V] = J 

I + 2V^) dx 
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(6.11) 

hiv] = J 

f ([Vf + IQV {V'f + dx , ... 

(6.12) 


It is an instructive exercise to compute these KdV integrals for a constant potential, and also for the 
Mathieu potential, to confirm the all-orders WKB analysis in Section 3. For example, for a constant 
potential, V = Vq, 


T /2 


2Vo3 




L Vie (4f;)3/2 (4f;)5/2 (4£;)7/2 


(6.13) 


which is the E —>■ +oo expansion of the exact result \/E — Vq. For the Mathieu system, to compare 
the Mathieu equation (2.1) with the Gelfand-Dikii form (6.7), we identify V O -^cosx, E u, 
and L -o- 2tt. The KdV integrals are simple to evaluate: 
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which leads to the expansion: 
Nh 


iv/6 +rr 1 1 11 

~Y ^ ““ 4(2u)3/2 “ 16 (2u)5/2 


- 


h? 


cosx 


64 (2u)7/2 


+ 15) - ... 


(6.14) 


(6.15) 


Note that the coefficients of the inverse powers of polynomials in of increasing order. This 

expansion (6.15) is a re-arrangement of the all-orders WKB expansion in (3.38). So, inverting (6.15) 
to express u as a function of h and a = Nh/2^ we recover the expansion (2.26) for u in the electric 
regime where a ^ 1. 

Turning now to the Lame system, Grosset and Veselov [84] considered the potential V = 2p,V, 
where V is the Weierstrass T^’-function (6.2), and ^ is a multiplicative scaling parameter to be specified 
below. Grosset and Veselov showed that the KdV integrals Ij+x\V] reduce to simple polynomials in pi, 
with coefficients expressed in terms of the Weierstrass invariants g 2 and gs, and ? 7 i/IK, where rji = C(If^) 
[84]: 


Ij+i[2gV] = Fj+i{p,) 


(6.16) 
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These polynomials, have interesting combinatorial properties, and were named the elliptic 

Faulhaber polynomials [84]^^. Thus, for this Weierstrassian form of the Lame equation, (6.7) with 
V{x) =V{x + iK';K,iK'), the high-energy Gelfand-Dikii expansion (6.8) can be written 


^ _ 1 \ " Pj+iilJ') 

2K "" ^ (4L;)J+V2 


E —>■ -boo 


(6.17) 


In order to facilitate the comparison with results in the physics literature concerning the Nekrasov 
prepotential in Af = 2* theories [30, 88, 89], we rewrite the coefficients of the elliptic Fauhaber 
polynomials in terms of the Eisenstein series 






(n,m)GZ2\(0,0) 


in -b rmY ’ 


(fc > 2). (6.18) 


where rjir) is the Dedekind eta function, and if 2 , E 4 and Eq can also be written as: 


E 2 =3 



m 

K 


E 4 


3 

4 



Ee 


27 

y 



(6.19) 


Thus, all coefficients of the elliptic Faulhaber polynomials are expressed in terms of just the first few 
Eisenstein series E2, E4 and Eq. In our conventions for the periods of the Weierstrass function, the 
modular parameter r is identified as 



( 6 . 20 ) 


For notational simplicity, we suppress the r (therefore k^) dependence of the Eisenstein series in the 
following equations. 

Then the first few elliptic Faulhaber polynomials F)+i(/i) are 

( 0 . 22 ) 

“ (s)°Y55(l-<256-9B2i54l(2M)’-12|E8-E2B4|(2f.)b (0.23) 

S^**'** “ (s)*iS (2rt* - 80 [B| - EjBe] (2 m)* + 96 [E| - EjBe] Plif) 

(6.24) 


Further properties of these elliptic Faulhaber polynomials F)_|_i(/r) are discussed in [84]. 

The Gelfand-Dikii expansion (6.17) can be viewed as a high-energy expansion of an all-orders 
Bohr-Sommerfeld relation, which can be inverted^^ to yield 

^^Not to be confused with the prepotential which we denote with the calligraphic letter J^. 

^^The name derives from a remarkable connection between the classical Faulhaber polynomials of Number Theory 
and the KdV integrals for soliton-like potentials [85], generalized to elliptic functions associated with periodic arrays of 
solitons [84]. In this sense, the polynomials in (6.14) could be referred to as trigonometric Faulhaber polynomials. 

^^It is interesting to note that the Lagrange inversion of series can be naturally expressed in terms of Young tableaux 
[22, 90]. 
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where the polynomials are simple combinations of the elliptic Faulhaber polynomials 


= -^[E 2 ]{ 2 t,) (6.26) 

(6.27) 

gsit^) = ^ ([2£;6 + ^E 2 E^ - 5Ei] { 2 fi)^ - 6 [Ee- E 2 E 4 ] (2^)^) (6.28) 

g^ifi) = ^ ([-35^2^ + 7EiEi + lOEj + 18E2Ee] {2^i)^ + 12 [-5F;| - 2E2Ee + lEjE^] (2/r)3 

uKjI 2 

+72 [-E 2 Ef, + El] (2^)2) (6.29) 


It is convenient to define the rescaled and shifted eigenvalue, absorbing the j = 0 term from the sum: 

2 


1 


2 K\^ E 2 , , 

- E+fi2,) 


(6.30) 


Then in terms of the action variable a = tiN/2, (6.25) becomes 


"~ 2 “ +2 


E 

1=1 ^ ^ 


1 2 , 




(6.31) 


It is a non-trivial check that with the scaling of (6.3) and (6.2), we identify 2/r e-i ^ 2 ^ 2 ; and in the 
1 0 limit we find that (6.31) does indeed reduce to the Mathieu large period expansion in (2.26). 
We can now compare the Bohr-Sommerfeld energy expansion (6.31) with the expansion of the 
Nekrasov-Shatashvili prepotential for the SU(2) Af = 2* theory [30, 86-89]^^: 


Ens ^ (®) 



no 

2TTi 


log(??(r)) 


1 ^ Hj 

Eki 2 ^ 2t+i j an 
j=i 


(6.32) 


where "Ho is expressed in terms of the scalar mass, m, and the Nekrasov deformation parameter in the 
Nekrasov-Shatashvili limit = ei as 


"Ho = 


4 


(6.33) 


The higher coefficients nj are polynomials in "Hq: 

ni(no) = l^mnl (6.34) 

^2 (Ho) = ^ ( [5 f;2" + Ei] nl - 3 [^4] nl n^) (6.35) 

nmo) = ([i75£;| + 84f;2£;4 + iieq] nt - 36 [7E2Ei + 8 Ee] nlh^ + i80[f;6] nlh'^) 

0U4oU 

(6.36) 

Using the differentiation properties of the Eisenstein series 

= , ni E,n-^EA-E,\ , +Te, = (6.37) 

^^Note that our normalization of T differs from the normalization in [30] as -Fhere = ~ 4 ^'^there* 
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and the identification 


2 ^ -O' 


Ho _ 


1 

4 


(6.38) 


we observe the remarkable fact that the expansion coefficients Hj are directly related to the polynomials 

Gj+i- 

, with 2/ro^ (6.39) 

We thus arrive at the Matone relation for the Af = 2* theory: 


■At=2* 


« - ^^'NS 



(6.40) 


The shift in u is just the shift in (6.30), required to have a smooth reduction to the Mathieu eigenvalue, 
as the Af = 2* theory reduces to the N = 2 theory in the m? —>■ oo limit, which is combined with 
—>■ 0 with finite.^® 

It is interesting to note that when the scalar mass m is related to the Nekrasov deformation 
parameter ei = h as 

m=(^j+^h => 2^ = j(j + l) (6.41) 

then if j is an integer the Lame system becomes finite gap [91], with the number of gaps equal to 
j. In this case, much more explicit descriptions of the spectral information can be given [84, 91], 
and there is also a precise strong/weak band/gap duality symmetry [81]. Notably for integer values 
of j, the associated density of states p{E)dE has an algebraic geometric meaning: it is an abelian 
differential of the second kind over a genus-j Riemann surface. Therefore the quantization condition 
ttN ^ f p{E)dE carries a clear geometric interpretation. 

Thus, in the large action spectral region of the Lame equation we identify the all-orders WKB 
Bohr-Sommerfeld expansion with the Gelfand-Dikii expansion, and the eigenvalue is identified with 
the scalar condensate moduli parameter of the N = 2* gauge theory. But, as in the Mathieu equation, 
this Bohr-Sommerfeld expression (6.31) for the energy is only part of the story. Identifying the action 
a with NH/2, the Bohr-Sommerfeld expression gives a (convergent) perturbative expression for the 
location of the gap high in the spectrum. However, physically it is clear that there are also non- 
perturbatively small gaps in the spectrum, and these are related to instantons. As in the Mathieu case 
discussed in Section 3, these gaps arise due to poles in the expansion coefficients: see the continued 
fraction expressions at http://dlmf.nist.gOv/29.3.iii. Since the structure is quite similar, we do 
not repeat all the steps here. Further discussion of the Lame system is deferred to a future publication. 


7 Conclusions 

In this paper we have applied resurgence and all-orders exact WKB to provide a complete description 
of the different spectral regions of the Mathieu and Lame systems, stressing their close connection with 
the low energy behavior oiAfi = 2 SUSY SU(2) gauge theories. Defining a’t Hooft parameter, X = N H, 
where N is the spectral level number, we associate the large A regime with the gauge electric regime. 

Effectively, we are subtracting the average value of the Weierstrass potential over one period, so that the first KdV 
integral, /i[V], vanishes, as it does for the Mathieu potential. 
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the small A regime with the gauge dyonic regime, and the A ^ 1 regime with the gauge magnetic 
regime. We have shown that exact WKB provides a complete description of all regimes, including 
all non-perturbative effects, and permits direct mappings between these sectors. Previous analyses 


based on all-orders Bohr-Sommerfeld relations described only part of the information, neglecting non- 


perturbative band or gap splittings. This analysis also shows that the relation between perturbative 
and non-perturbative contributions is radically different in the different regimes. The familiar rela¬ 
tion between divergent perturbative expansions and non-perturbative effects must be generalized to 
accommodate the large A regime where perturbative expansions are convergent, and non-perturbative 
effects are associated with poles of expansion coefficients. This generalizes the resurgent trans-series 
analysis of [9, 34] to large A. We provide a physical analogy of this change of non-perturbative behav¬ 
ior, in the transition from tunneling pair production to multi-photon pair production in the Schwinger 
effect in a time-dependent electric field, as the adiabaticity of the field changes. This physics is natu¬ 
rally described by worldline instantons. Since the Nekrasov deformation parameters have the physical 
interpretation of constant graviphoton fields, the SUSY gauge theory significance of these worldline 
instantons, and their associated non-perturbative effects, should be further understood. 

The reinterpretation of the spectral problem in the language of the Nekrasov partition function 
leads naturally to a simple proof of a quantum mechanical resurgence relation [9] that shows that 
all non-perturbative information of the trans-series eigenvalues is subtly encoded within perturbation 
theory, which is an extreme form of resurgence behavior. We also give a geometrical interpretation of 
this fact in terms of the Riemann bilinear identity and the Picard-Fuchs equation, complementary to 
work on quantum geometry [25, 30, 32] and Whitham dynamics [33]. We demonstrate a direct relation 
between the all-orders Bohr-Sommerfeld relation and the Gelfand-Dikii expansion, based on the KdV 
invariants, and show explicitly how the JV = 2* theory reduces to the N = 2 system. Future work will 
develop these ideas further [77]. 

Acknowledgements. We thank C. Bender, O. Costin, A. Gorsky, M. Marino, N. Nekrasov, M. 
Unsal and A. Vainshtein for discussions and correspondence. GB thanks the DOE for support through 
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oretical Physics Institute at Friedrich-Schiller Universitat, Jena, where much of this work was done. 
GD also thanks the DFG for a Mercator Professor Grant, and the US-German Fulbright Gommission 
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A Simple Analog of Uniform Asymptotic Behavior 

As discussed in [8, 92], certain features of the divergence of the perturbative expansion for the ground 
state energy of the Mathieu system are captured by the zero-dimensional partition function 



(A.l) 


The relation to the quantum mechanical spectral problem arises because this zero-dimensional partition 
function gives the resummed leading derivative expansion contribution to the heat kernel expansion 
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tre“^*, from which the resolvent and spectrum can be extracted. The perturbative expansion in 
(A.l) is asymptotic, with non-alternating factorial large order behavior of the expansion coefficients, 
Cn ^ {n— l)!/ 7 r. The series is non-Borel-summable and should instead be represented by a trans-series, 
which has just two exponential terms since Z{g) satisfies a second order differential equation [2]:^® 


Z{9) 




69 


r(r 


n =0 


TT 2"'n! 


'± 


OO 


^ 2 ^ 


X 2 ) g" 

^ 7r2"n! ^ 

n—Q 


(A.2) 


where — f -|- d < Tphg < | 7 r — d. This trans-series expression contains all information about the 
function, including the Stokes phenomenon. In particular, the second sub-dominant part is required 
in order to be consistent with the analytic continuation connection formula for the Bessel functions: 


A:o(e^*’"z) = Kq{z) =F *7r/o(z) 


(A.3) 


Note that in the large g limit this function has a convergent “strong-coupling” expansion: 


^ 1 1 

Z{g) > -9 ) 9 +00 

The dependence on the level number N can be modeled by considering instead 


(A.4) 


Zn{9) = - 

\9 


\/27r/5 


e'n'^Cn{N) g'^ , 5 ^ 0 + (TV fixed) 


(A.5) 


n—0 


where the series is again divergent and non-Borel-summable, with c„(N) ^ (_ 1 )-^ > sf large 

perturbative order n, as in (A.l). But we can also consider the large N limit with g fixed: 


Zn{9) 


V2 ttN \2 Ng 


N 


, N —>• -boo {g fixed) 


(A. 6 ) 


which has the same form as the gap splitting (2.30) high in the spectrum where N 3> 1/g. On the 
other hand, there is also a uniform expansion where N ^ 00 and 5 —> 0 such that Ng is kept fixed. 
The uniform expansion valid for all values of the’t Hooft parameter X = N g is: 
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1 exp 
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Ng 


N 




(«*+ X) 


1/4 


\f^ 


1 + a/1 + 




(Ngy 

N —>■ -boo (0 < A^ 5 < 00 ) 


n =0 


i-b 


(Ng)^ 


(A.7) 


where Un is a (known) polynomial of degree 3n [http://dlmf.nist.gov/10.41.E3]. This uniform 
expression interpolates smoothly between the two extreme limits and describes the function in the 
intermediate region, the analog of the “magnetic region”. 

^®For the Lame potential, Z(g) satisfies a third-order differential equation, and this fact is reflected in the appearance 
of three different exponential terms, which can be identified with the perturbative vacuum and both real and complex 
(’ghost’) instantons [35]. 
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B WKB coefficients 


In this appendix we list the first five set of coefficients k 


(n) 
k ’ 


that enter into the relation 


n 

an{u) = 


(n) fc d"+'^«0(M) 

k du'^+^ 


(B.l) 


k=0 

that connects the higher order WKB cycles to the leading order one. n denotes the coefficient of 
in the WKB expansion. The method for calculating them is explained in Section 3.4. 


n 4, /fc — >■ 

0 

1 

2 

3 

4 

5 

1 

1/48 

1/24 





2 

5/1536 

1/192 

7/5760 




3 

41/57344 

153/143360 

79/215040 

31/967680 



A 

15229 

9539 

517 

13 

127 



70778880 

30965760 

4128768 

716800 

154828800 


5 

484249 

5049503 

8430053 

780341 

61729 

73 

5813305344 

43599790080 

163499212800 

81749606400 

81749606400 

3503554560 
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